Critical currents in superconductors with quasiperiodic pinning arrays: 
One-dimensional chains and two-dimensional Penrose lattices 
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We study the critical depinning current Jc, as a function of the applied magnetic flux for 
quasiperiodic (QP) pinning arrays, including one-dimensional (ID) chains and two-dimensional (2D) 
arrays of pinning centers placed on the nodes of a five-fold Penrose lattice. In ID QP chains 
of pinning sites, the peaks in Jc($) are shown to be determined by a sequence of harmonics of 
long and short periods of the chain. This sequence includes as a subset the sequence of successive 
Fibonacci numbers. We also analyze the evolution of Jc('&) while a continuous transition occurs from 
a periodic lattice of pinning centers to a QP one; the continuous transition is achieved by varying 
the ratio 7 = ag/a^ of lengths of the short as and the long segments, starting from 7 = 1 for a 
periodic sequence. We find that the peaks related to the Fibonacci sequence are most pronounced 
when 7 is equal to the "golden mean". The critical current Jc($) in QP lattice has a remarkable 
self-similarity. This effect is demonstrated both in real space and in reciprocal fc-space. In 2D 
QP pinning arrays (e.g., Penrose lattices), the pinning of vortices is related to matching conditions 
between the vortex lattice and the QP lattice of pinning centers. Although more subtle to analyze 
than in ID pinning chains, the structure in Jc{^) is determined by the presence of two different kinds 
of elements forming the 2D QP lattice. Indeed, we predict analytically and numerically the main 
features of Jc($) for Penrose lattices. Comparing the Jc's for QP (Penrose), periodic (triangular) 
and random arrays of pinning sites, we have found that the QP lattice provides an unusually broad 
critical current Jc{^), that could be useful for practical applications demanding high Jc's over a 
wide range of fields. 

PACS numbers; 74.25.Qt 
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I. INTRODUCTION 

Recent progress in the fabrication of nanostructures 
has provided a wide variety of well-controlled vortex- 
confinement topologies, including different types of regu- 
lar pinning arrays. A main fundamental question in this 
field is how to drastically increase vortex pinning, and 
thus the critical current Jc, using artificially-produced 
periodic arrays of pinning sites (APS). These periodic 
APS have been extensively used for studying vortex pin- 
ning and vortex dynamics. In particular, enhanced Jc 
and commensurability effects have been demonstrated in 
superconducting thin films with square and triangular 
arrays of sub-/^m holes (i.e., antidots) H E H i, H H . 
Moreover, blind antidots (i.e., holes which partially per- 
forate the film to a certain depth) or pinning ar- 
rays with field-dependent pinning strength ||8|, provide 
more flexibility for controlling properties such as pinning 
strength, anisotropy, etc. The increase and, more gener- 
ally, control of the critical current Jc in superconductors 
by its patterning (perforation) can be of practical im- 
portance for applications in micro- and nanoelectronic 
devices. 

A peak in the critical current Jc($), for a given value 
of the magnetic flux, say $1, can be engineered using 
a superconducting sample with a periodic APS with a 
matching field Hi = ^i/A (where A is the area of the 
pinning cell), corresponding to one trapped vortex per 



pinning site. However, this peak in Jc($), while useful 
to obtain, decreases very quickly for fluxes away from <f>i. 
Thus, the desired peak in Jc($) is too narrow and not 
very robust against changes in <i>. It would be greatly 
desirable to have samples with APS with many periods. 
This multiple-period APS sample would provide either 
very many peaks or an extremely broad peak in Jc($), 
as opposed to just one (narrow) main peak (and its har- 
monics). We achieve this goal (a very broad Jc{^)) here 
by studying samples with many built-in periods. 

The development of new fabrication technologies for 
pinning arrays with controllable parameters allows to 
fabricate not only periodic (square or triangular) but also 
more complicated quasiperiodic (QP) arrays of pinning 
sites, including Penrose lattices [9l[l0l[lll|. 

The investigation of physical properties of QP systems 
has attracted considerable interest including issues such 
as band structure and localization of electronic states in 
two-dimensional (2D) Penrose lattice 12, 13], electronic 
and acoustic properties of one-dimensional (ID) QP lat- 
tices 0, 0, superconducting-to-normal phase bound- 
aries of 2D QP micronetworks 0, IH El IH , QP 
semiconductor heterostructures and optical superlattices 
|2ll |. soliton pinning by long-range order a-^d pulse 
propagation l23!| in QP systems. Moreover, increasing 
and, more generally, controlling the critical current in 
superconductors by its patterning (perforation) can be 
of practical importance for applications in micro- and 
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nanoelectronic devices. 

The original tiling has been studied in Ref. (21] ■ The in- 
flation (or production) rules of "finite Penrose patterns" 
generated by repeated application of deflation and rescal- 
ing have been found, which show a definite hierarchical 
structure of the Penrose patterns . 

The electronic and acoustic properties of a one- 
dimensional quasicrystal have been studied in Refs. ll^ 
IT^ . It has been shown, in particular, that there exist two 
types of the wave functions, self-similar (fractal) and non- 
self-similar (chaotic) which show "critical" or "exotic" 
behavior ■ By both numerical (non-perturbative) and 
analytical (perturbative) approaches, it has been demon- 
strated 14, 15] that the phonon and electronic spectra of 
ID quasicrystals exhibit a self-similar hierarchy of gaps 
and locahzed states in the gaps. The existence of gaps, 
and gap states, in QP GaAs-AlGaAs superlattices has 
been predicted and found experimentally. 

Along with studying the structural, electronic and 
acoustic properties of QP structures, considerable 
progress has been reached in understanding the su- 
perconducting properties of 2D quasicrystalline arrays 
111 E E lia IIS- The effect of frustration, induced 
by a magnetic field, on the superconducting diamagnetic 
properties has been revealed and the superconducting- 
to-normal phase boundaries, Tc{H), have been calcu- 
lated for several geometries with quasicrystalline order, 
in a good agreement with experimentally measured phase 
boundaries [TtI Il9j ] . A comprehensive analysis of super- 
conducting wire networks including quasicrystalline ge- 
ometries and Josephson-junction arrays in a magnetic 
field has been presented in Refs Ej 113 • An analytical 
approach |l9l l2Cljl was introduced to analyze the struc- 
tures which are present in phase diagrams for a number 
of geometries. It has been shown that the gross structure 
is determined by the statistical distributions of the cell 
areas, and that the fine structures are determined by cor- 
relations among neighboring cells in the lattices. The ef- 
fect of thermal fluctuations on the structure of the phase 
diagram has been studied [l^ by a cluster mean-fleld 
calculation and using real-space renormalization group. 

In this paper, we study another phenomenon related 
to superconducting properties of quasiperiodic systems, 
namely, vortex pinning by ID QP chains and by 2D ar- 
rays of pinning sites located at the nodes of QP lattices 
(e.g., a flve-fold Penrose lattice). It should be noted that 
in superconducting networks the areas of the network pla- 
quettes play a dominant role However, 
for vortex pinning by QP pinning arrays, the specific ge- 
ometry of the elements which form the QP lattice and 
their arrangement (and not just the areas) are impor- 
tant, making the problem complicated. 

In Sec. II we introduce the model used for describing 
vortex dynamics in QP pinning arrays and for determin- 
ing the critical depinning current, Jc, which is analyzed 
for different quasicrystalline geometries. 

The pinning of vortices by a ID QP chain of pinning 
sites (i.e., Fibonacci sequence) is discussed in Sec. III. We 



consider a continuous transition from a periodic to the 
QP chain of pinning sites and we monitor the correspond- 
ing changes in the critical current Jc as a function of the 
applied magnetic flux <&. A remarkable self- similarity of 
Jci^) is demonstrated in both real space and in recipro- 
cal fc-space. 

Section IV studies the pinning of flux lattices by 2D 
QP pinning arrays including the flve-fold Penrose lat- 
tice. We analyze changes of the function Jc{^) during 
a continuous transition from a periodic triangular lat- 
tice of pinning sites to the Penrose lattice. Based on 
detailed considerations of the structure and specific local 
rules of construction of the Penrose lattice, we predict 
the main features of the function Jc{^)- Numerical sim- 
ulations with different finite-size Penrose lattices confirm 
the predicted main features for large-size lattices, that is 
important for possible experimental observations of the 
revealed quasiperiodic features. We also obtain analyt- 
ical results supporting our conclusions. Moreover, we 
also discuss the changes in the critical current by adding 
either a "quasiperiodic" modulation or random displace- 
ments to initially periodic pinning arrays. 

II. MODEL 

We model a three-dimensional (3D) slab infinitely long 
in the z-direction, by a two-dimensional (2D) (in the xy- 
plane) simulation cell with periodic boundary conditions, 
assuming the vortex lines are parallel to the cell edges. 
To study the dynamics of moving vortices driven by a 
Lorentz force, interacting with each other and with pin- 
ning centers, we perform simulated annealing simulations 
by numerically integra ting the overdamped equations of 
motion (see, e.g., Ref. IM IH HI 113 ) : 

77V, = f , = r+fr+ff +f^ (i) 

Here, fi is the total force per unit length acting on vortex 
i, f™ and f"^ are the forces due to vortex- vortex and 
vortex-pin interactions, respectively, is the thermal 
stochastic force, and ff is the driving force acting on the 
i-th vortex; rj is the viscosity, which is set to unity. The 
force due to the interaction of the i-th vortex with other 
vortices is 

r = E fo K,0-^^^) , (2) 

where Ny is the number of vortices, Ki is a modified 
Bessel function, A is the magnetic field penetration depth, 
= (rj - Tj)/ \ri-rj |, and 

f ^0 
■"^ 87r2A3 ■ 

Here $o — hc/2e is the magnetic flux quantum. It is 
convenient, following the notation used in Ref. [23 . l25l 
m 113] , to express all the lengths in units of A and all 
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the fields in units of ^q/}?. The Bessel function Ki{r) 
decays exponentially for r greater than A, therefore it is 
safe to cut off the (negligible) force for distances greater 
than 5A. The logarithmic divergence of the vortex- vortex 
interaction forces for r — > is eliminated by using a cutoff 
for distances less than O.IA. 

Vortex pinning is modeled by short-range parabolic po- 
tential wells located at positions r^-^' . The pinning force 
is 
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where Np is the number of pinning sites, fp is the maxi- 
mum pinning force of each potential well, rp is the range 
of the pinning potential, Q is the Heaviside step function. 
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The temperature contribution to Eq. Q is represented 
by a stochastic term obeying the following conditions: 



(4) 



and 



Throughout this work, we use narrow potential wells as 
pinning sites, characterized by rp = 0.04A to O.IA. Our 
calculations show that, for the parameters used, the de- 
pendence of the critical current Jc{^) defined accord- 
ing to Eq. lO, is in good agreement to that based on 
the above general definition of Jc (which involves an ad- 
justable parameter Jmin)- The advantages of using Jc 
defined by Eq. © are the following: it (i) does not in- 
volve any arbitrary threshold Jmin and (ii) is less CPU- 
time consuming, allowing the study of very large-size lat- 
tices. Moreover, the goal of this study is to reveal spe- 
cific matching effects between a (periodic) vortex lattice 
and arrays of QP pinning sites, and to study how the 
quasiperiodicity manifests itself in experimentally mea- 
surable quantities (Jc, Ni^^^/N^) related to the vortex 
pinning by QP (e.g., the Penrose lattice) pinning arrays. 



III. PINNING OF VORTICES BY A ID 
QUASIPERIODIC CHAIN OF PINNING SITES 

In this section we study the pinning of vortices by one 
dimensional (ID) QP chains of pinning sites. 



{fl{t)ff{t'))=2f^kBT6,,6{t~t'). 



(5) 



The ground state of a system of moving vortices is ob- 
tained as follows. First, we set a high value for the tem- 
perature, to let vortices move randomly. Then, the tem- 
perature is gradually decreased down to T = 0. When 
cooling down, vortices interacting with each other and 
with the pinning sites adjust themselves to minimize the 
energy, simulating the field-cooled experiments [s^l . 

In order to find the critical depinning current, Jc, we 
apply an external driving force gradually increasing from 
/d = up to a certain value fd = /J, at which all the 
vortices become depinned and start to freely move. For 
values of the driving force just above /J, the total cur- 
rent of moving vortices J ~ (v) becomes nonzero. Here, 
(v) is the normalized (per vortex) average velocity of all 
the vortices moving in the direction of the applied driving 
force. In numerical simulations, this means that, in prac- 
tice, one should define some threshold value Jmin larger 
than the noise level. Values larger than Jmin are then con- 
sidered as nonzero currents. However, instead using this 
criterion-sensitive scheme, we can use an alternative ap- 
proach based on potential energy considerations. In case 
of deep short-range potential wells, the energy required 
to depin vortices trapped by pinning sites is proportional 
to the number of pinned vortices, Ni^\ Therefore, in 
this approximation we can define the "critical current" 
as follows: 



Jc($) = Jo 



(6) 



where Jq is a constant, and study the dimensionless value 
J^ = Jc/Jq (further on, the primes will be omitted). 



A. ID quasiperiodic chain 

As an example of a ID QP chain, or ID quasicrystal, 
a Fibonacci sequence is considered, which can be con- 
structed following a simple procedure: let us consider 
two line segments, long and short, denoted, respectively, 
by L and S. If we place them one by one, we obtain an 
infinite periodic sequence: 



LSLSLSLSLSLSLSLS . 



(7) 



A unit cell of this sequence consists of two elements, L 
and S. In order to obtain a QP sequence, these elements 
are transformed according to Fibonacci rules as follows: 
L is replaced by LS, S is replaced by L: 



As a result, we obtain a new sequence: 
LSLLSLLSLLSLLSL . 



(8) 



(9) 



Iteratively applying the rule (|SJ to the sequence lO, 
we obtain, in the next iteration, a sequence with a five- 
element unit cell (LSLLS), then with an eight-element 
unit cell {LSLLSLSL), and so on, to infinity. For the 
sequence with an n-element cell, where n — s- cx), the ratio 
of numbers of long to short elements is the golden mean 
value, 



(1 + \/5)/2. 



(10) 



The position of the nth point where a new element, either 
L or S, begins is determined by pTl |: 



(11) 
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FIG. 1: (a, b) Dimensionless critical depinning current Jc, 
as a function of the applied magnetic flux <1>, in a ID periodic 
chain of pinning sites, for a cell containing 25 pinning sites, 
Np = 25 and fp/fo = 2.0. The indicated fluxes $i, <1>2, 
and $3 correspond to the first, second, and third matching 
fields. The function Jc($) is shown for two different values 
of the pinning site radius, rp — 0.1 A (a), and rp — 0.04A 
(b). When Vp, decreases, the main commensurability peaks 
become sharper, as shown in (a) and (b). 



B. Pinning of vortices by a ID periodic chain of 
pinning sites 

We start with a periodic chain of pinning sites, which 
can be considered, in the framework of the above scheme, 
as a limiting case of a "QP" chain with 7 = 1, i.e. as ~ 
flL = 1. In Fig. la, the critical current Jc N'f'' /Ny 
is shown as a function of the applied magnetic flux <I> 
for /p//o = 2 and for a pinning site radius = O.IA. 
Sharp peaks of the function Jc($) correspond to match- 
ing fields. Since the dimensionless critical current we 
plot is that per vortex and is proportional to the num- 
ber of pinned vortices and inversely proportional to the 
total number of vortices (Eq ©), therefore, the mag- 
nitude of the peaks versus ^ ^ Ny decreases as 1/$. 
Then the maximum heights of the peaks are: Jc($i) = 1, 
Jc(2$i) = 0.5, and Jc(3$i) = 0.33. Note that these val- 
ues are obtained provided each pinning site can trap only 
one vortex, which is justified for the chosen radius of the 
pinning site and for the vortex densities considered. In 
addition, there are weak wide maxima corresponding to 
$i/2 and other "subharmonics" , i.e. <i>i-|-$i/2 (Fig. la), 
where i is an integer. For very small values of $i, the 
vortex density is very low, and vortices almost do not 
interact with each other. As a result, in the ground state 
they all become trapped by pinning sites, and Jc is max- 
imal for small $i. 

For smaller radii of the pinning sites, rp = 0.04A 
(Fig. lb), the Jc{^) peaks corresponding to matching 
fields become sharper because for smaller values of rp, 
it is more difficult to fulfill the commensurability con- 
ditions. Any features around the main peaks, including 
subharmonics, are suppressed. Note also a "parity effect" 
takes place in this case: since the number of pinning sites 
per cell is odd {Np = 25 in Fig. 1), therefore Jc{^) peaks 
are suppressed for <i>i/2 and for even values of i in the 
sequence $i + $1/2. 



where [x\ denotes the maximum integer less or equal to 
X. Equation (|llll corresponds to the case when the Fi- 
bonacci sequence has a ratio 7 = as/a^ of the length of 
the short segment, as, to the length of the long segment, 
a^. Ratios 7 other than 1/r correspond to other chains 
which are all QP. Along with the golden mean value of 
7 — l/r, we use in our simulations 7's varying in the 
interval between and 1: < 7 < 1, when analyzing a 
continuous transition from a periodic to a QP (Fibonacci 
sequence) pinning array. 



To study the critical depinning current Jc in ID QP 
pinning chains, we place pinning sites on the points where 
L OT S elements of the QP sequence link to each other. 
Therefore, the coordinates of the centers of the pinning 
sites are defined by Eq. (fTT)) with j = as/aL- 



IV. PINNING OF VORTICES BY A ID 
QUASIPERIODIC CHAIN OF PINNING SITES: 
GRADUAL EVOLUTION FROM A PERIODIC 
TO A QUASIPERIODIC CHAIN 

Let us consider a QP chain of pinning sites with spac- 
ings between pinning sites given by ol (long) and as 
(short). The long and short segments alternate accord- 
ing to the Fibonacci rules forming a Fibonacci sequence. 
The number of pinning sites per cell coincides with the 
number of elements of this sequence per unit cell. It is 
natural to take a chain with a number of elements equal 
to one of the successive Fibonacci numbers as a ID cell, 
although in principle it could be of any length. We im- 
pose periodic boundary conditions at the ends of the cell. 

The larger cell we take, the closer we are to describ- 
ing a truly QP structure. However, it turns out, that 
even a finite part of a QP system (ID chain or 2D QP 
lattice) provides us with reliable information concerning 
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FIG. 2: (Color) Critical depinning current Jc, as a func- 
tion of the applied magnetic flux (for convenience, shown as 
a function of the number of vortices, Nv ~ $), in a ID QP 
chain of pinning sites with fp/fo — 1, rp — 0.1 A, for chains 
characterized by different ratios of the lengths of short as to 
long aL spacings, 7 — as/aL- For small deviations of the 
chain from a periodic chain (7 = 0.9 to 1), commensurability 
peaks are similar to those shown in Fig. 1. For intermedi- 
ate values of 7 (7 = 0.2 to 0.8), peaks are determined by 
a sequence of harmonics of numbers of long and short peri- 
ods of the chain, which includes the sequence of successive 
Fibonacci numbers, most pronounced for 7 = 1/r « 0.618, 
where t — {1 + \/5)/2 ~ 1.618 is the golden mean. For very 
small 7 (e.g., 7 = 0.1), the QP chain effectively becomes peri- 
odic but with the number of pinning sites equal to the number 
of long periods in the chain. 

properties of the whole system. This is based on the 
structural self- similarity of QP systems. These prop- 
erties are studied here for the critical depinning current 
Jc and have also been demonstrated for other physical 
phenomena [ll|ll[llll8, 19, 20, 28]. 

Figure 2 shows the evolution of Jc as a function of the 
number of vortices, ~ 4>, for various values of the 
parameter 7. The top curve represents the limiting case 
of a periodic chain (7 = 1.0) with typical peak structure 
discussed above (Fig. 1). The chain contains Np = 21 
pinning sites. As a result, commensurability peaks ap- 
pear at Ny = 21, 42, 63, etc., i.e. multiples of 21. A small 
QP distortion of the chain does not appreciably affect the 
peak structure of the function Jc{Ny) (7 = 0.9). When 
the deviation of the factor 7 from unity becomes larger 
(7 = 0.8, 0.7), commensurability peaks for Ny — 42, 63, 
and other multiples of 21 decrease in magnitude. At the 
same time, new peaks appear at Ny = 55, 76 (= 55-1-21), 
97 (= 76 + 21) for 7 = 0.8. Then, with further decrease 
of 7 (7 = 0.7), these peaks remain {Ny = 55 even grows 
in magnitude), and a new peak at iV„ = 34 arises. 

For the golden mean-related value of 7 = 1/r « 0.618, 
we obtain a set of peaks, which are "harmonics" of the 
numbers of long and short periods of the chain (or recip- 
rocal lengths ul and as), i.e. 

4aks.. = A,N, + B,Ns = ^ + ^, (12) 

as aL 

where Nl and Ns are the numbers of long and short el- 
ements, respectively; Ai (A^) and Bi [B'^) are generally 
(positive or negative) multiples or divisors of N^ and 
Ns\ the upper index "QP" denotes "quasiperiodic" . It is 
easy to see that this set includes as a subset the sequence 
of successive Fibonacci numbers. In particular, the fol- 
lowing well-resolved peaks of the function Jc{Ny) appear 
for 7 = 1/t « 0.618: Ny = 13 (= Nl = 13, 13 is a Fi- 
bonacci number (FN)); Ny = 17 (= {2Nl+Ns)/2, where 
{2Nl + Ns) = 34 is a FN); Ny = 21 (= Np ^ Nl + NS, 
21 is a FN); Ny w 27-28 (= {^l + 2iVs)/2, where 
{3Nl + 2Ns) = 55 is a FN); Ny = 34 {2Nl + Ns); 
Ny « 44-45 (= (5iVL-l-3A^s)/2, where {5Nl + SNs) = 89 
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FIG. 3: (Color) (a) The critical depinning current Jc, as a 
function of the number of vortices, ~ $, for different ID 
QP chains, Np = 21 (red bottom line), Np — 34 (blue line), 
Np — 55 (green line), and Np = 89 (dark blue top line), and 
the same 7 = as/aL = 1/t. The parameters used here are: 
fp/ fo = 1.0 and rp — O.IA. Independently of the length of 
the chain, the peaks for all of the curves include as a subset 
the sequence of successive Fibonacci numbers (indicated by 
blue arrows in the horizontal axis) and their sub-harmonics, 
(b) Jc{Nv) for a long chain Np = 144 and the same 7 = 
1/r. Notice that now A'^ is much larger than in (a), (c) The 
function Jc ($/$!) for the same set of ID chains (using the 
same colors), normalized by the number of pinning sites for 
each chain. The curves for chains with different A^p's display 
the same sets of peaks, namely, at $/$i = 1 (first matching 
field) and "l>/"l>i — 0.5, as well as at the golden-mean-related 
values: $/$i = r/2, $/$i = {r + l)/2 = t^/2, $/$i = r, 
^/^i = (r+r2)/2 = = = t+1, <I>/$i = + 

Therefore, similar sets of peaks are obtained for both cases: 
for the curves plotted in the same scale (a), (b); for the curves 
plotted in individual scales, i.e., normalized on the number of 
pinning sites in the chain (c) . This behavior demonstrates the 
self-similarity of Jc{^)- 



is a FN); iV,„ = 55 (= {3Nl + 2Ns), 55 is a FN); iV„ = 68 
(= (47Vi + 2Ns)); = 89 (= {5Nl + 3Ns), 89 is a 
FN). In summary, the most pronounced peaks are at: 
Ny = 13, 17, 21, 34, 55, 89, which (except the point 
Ny = 17 = 34/2) form a sequence of successive FNs. 

The above QP peaks only slightly degrade at 7 = 0.6. 
However, when the length of the long segment be- 
comes twice the length of the short one 05, i.e. 7 — 0.5, 
sharp commensurability peaks appear which are related 
to the small segment of the chain with length as — 0,^/2. 
Namely, we obtain peaks at A'^i, = 34 ~ {2Nl + Ns) 
and at other values of Ny which are sub-harmonics of 
2Nl + Ns: Ny = 17, 34, 68. Other peaks, in particu- 
lar, those related to the Fibonacci sequence, are much 
less pronounced for 7 — 0.5. For 7 — 0.4, the peaks 
are at: Ny = 13, 17, 34, 47 (= 3Nl + Ns), 68, and 81. 
A very strong (recall that the maximum amplitude of 
the peak is ~ 1/Ny) peak at Ny = 81 is a "resonance" 
peak corresponding to the ratio of Bi/Ai = 2/5 = 7, i-^- 
57Vl + 2A^s = 81. When 7 = 0.3, the "resonance" ratio 
Bi/Ai ~ 0.3, therefore, a strong peak appears for the 
"nearest to 0.3" value B,/Ai = 0.33 : ZNl + Ns ^ 47, 
and also for 6Nl + 2Ns = 94. Another "close to 0.3" 
value is Bi/Ai — 0.25, which is responsible for the peaks 
at Ny = 30 {2Nl+Ns/2) and Ny = 60 {ANl+Ns). Also, 
peaks at Ny = 13, 17 and 55 are present. The "resonant" 
peak for 7 = 0.2 (i.e., for the ratio 7 = Bi/Ai = 0.2), 
appears at Ny = 73 (= 5iVi -f Ns)- The closest neigh- 
boring peaks are at Ny = 60 (= ANl + Ns) and Ny = 86 
(= QNl + Ns), [also: Ny = 43 (= 3iVL + iVs/2)], charac- 
terized by ratios Bi/Ai = 0.25 and Bi/Ai — 0.17, corre- 
spondingly. Finally, for 7 = 0.1 we arrive at the situation 
when we have an almost periodic chain but with a period 
a different from that for 7 = 1.0: a' « since as <IC aL- 
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The number of pinning sites becomes Np — Nl — 13, 
and we obtain commensurability peaks at the positions 
which are multiples of Np = 13, i.e., at Ny = 13, 26, 39, 
52, which is typical for periodic chains of pinning sites. 

In order to demonstrate that the above analysis is gen- 
eral and reveals the QP features independently of the 
length of a specific chain of pinning sites, let us compare 
results from the calculation of Jc{Ny) for different chains. 
In Fig. 3a, the function Jc{Ny) is shown for four different 
ID QP chains, Np = 21, Np = 34, Np = 55, and Np = 89, 
and the same 7 = 1/r. Figure 3a clearly shows that the 
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positions of the main peaks in Jc, i.e. those corresponding 
to a Fibonacci sequence, and other peaks whose positions 
are described by Eq. I|12(l . to a significant extent, do not 
depend on the length of the chain. The peaks shown in 
Fig. 3 form a Fibonacci sequence: A^^, = 13, 21, 34, 55, 89, 
144, and other "harmonics": 7V„ = 17, 27-28 (= 55/2), 
44-45 {— 89/2), etc. At the same time, longer chains al- 
low to better reveal peaks for larger Fibonacci numbers. 
Thus, for chains with Np = 144 (Fig. 3b) peaks at the 
next Fibonacci numbers are pronounced: Ny = 144, 233, 
377. 

While the curves for different chains are plotted in 
Fig. 3a in the same scale. Fig. 3c shows these curves 
in individual scales. Namely, we rescale each Jc by nor- 
malizing each Jc by the number of pinning sites in each 
curve. Thus, $1 corresponds to iV„ = 21 for the chain 
with Np = 21, to Ny = 34 for the chain with Np = 34, 
etc. After this rescaling, the Jc curves approximately fol- 
low each other and have pronounced peaks for the goldcn- 
mean-related values of <i>/$i, as shown in Fig. 3c. For 
example, $/<I>i — t corresponds to Ny = 34 for the chain 
with Np = 21, to Ny — 55 for the chain with Np = 34, to 
Ny = 89 for Np = 55, to Ny = 144 for Np = 89, and to 
Ny = 233 for the chain with Np = 144. Note that these 
peaks (i.e., corresponding to the golden mean) are most 
pronoimced for each chain in Fig. 3a, b. 

Therefore, the same peaks of the function Jc($) for 
different chains are revealed before and after rescaling. 
This means that the function Jc($) for the ID QP chain 
is self- similar. Below, we demonstrate the revealed self- 
similarity effect in a reciprocal /c-space. 



A. Fourier-transform of the vortex distribution 
function on a ID periodic chain of pinning sites: 
Self-similarity effect 

As we established above, the revealed QP features 
(e.g., peaks of the function Jc(<i> ^ Ny) ) (i) are inde- 
pendent of the length of the chain and (ii) the longer 
chain we take, the more details (e.g., "subharmonics" ) of 
QP features can be observed. 

This result is related to an important property of QP 
systems, self-similarity, which could be better understood 
by analyzing the Fourier-transform of the distribution 
function of the system of vortices pinned on a QP array. 

In /c-space, the distribution function Fy^"^ (q) of a sys- 
tem of Ny pinned vortices can be represented as the 

inverse Fourier-transform of the ID distribution function 

Ap) 



Fy (n) of the vortices in real space: 



J2 F^i^Hn) exp{-27rign/Ar,(P)}. (13) 




Figure 4 shows the Fourier-transform of a system of 
pinned vortices for a 7 equal to the inverse golden mean: 
1/7 = ttL/as = T = {1 + \/5)/2. The plots shown in 



FIG. 4: The self-similar Fourier-transform of the distribution 
function (vortex density) of the system of — 144 vortices 
pinned on a QP array, for 7 = 1/r. The portion limited by 
the two arrows in (a) is successively magnified several times 
and the corresponding results shown in (b), (c) and (d). 



Figs. 4a to 4d are obtained according to the following 
rule. The portion of the horizontal axis which is limited 
by the two arrows in Fig. 4a is rescaled and shown in 
Fig. 4b. In the same way, the portion limited by the 
two arrows in Fig. 4b is rescaled and shown in Fig. 4c. 
Fig. 4d is obtained following the same procedure. Note 
that each subsequent scaling is accompanied by flipping 
the direction of the g-axis to the opposite direction. A 
similar property is clear from the experimental diffraction 
patterns of quasicrystals fllj . The pentagons of Bragg 
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FIG. 5; Fourier-transform of the distribution function of the 
vortices {Nv = 144) interacting with a QP array of pinning 
sites, for different values of 7 = 1.0 (a), 0.8 (b), 0.5 (c), 0.1 (d). 
The value 7 = 1 corresponds to the limit of a periodic chain. 
Varying 7 from one, we introduce quasiperiodicity (the most 
pronounced for 7 = 1/r « 0.618) in the chain. The 7 = 1 case 
recovers the periodic limit (with another period). The main 
central sharp peak, corresponding to the periodic chain used 
for (a), continuously transforms — through the set of self- 
similar patterns ((b) and (c)) corresponding to QP chains 
(see Fig. 4) — to another peak (d) produced by a periodic 
chain with the number of sites equal to the number of long 
segments of the initial periodic chain with 7 = 1.0 shown in 
(a). 



peaks have smaller pentagons inside them, which are in- 
verted. As seen in Fig. 4, each subsequent subdivision 
leads to a subset of peaks similar to the entire set of 
peaks. 

This analysis clearly demonstrates the self-similarity of 
the distribution function of the vortices pinned on a QP 
ID array of pinning sites. Similarly to the observed be- 



havior of the function Jc{^) when increasing the length 
of the chain of pinning sites, the Fourier-transform of the 
distribution function of the vortices pinned on a QP ar- 
ray reproduce its main features (peaks) in a self-similar 
way, when increasing the range in fc-space, and simulta- 
neously acquires a more elaborate structure with smaller 
self-similar peaks. 

As we discussed above, the main commensurability 
peaks evolve from a perfectly periodic set of the type 
<I>i = m<i>i, where m > 1 is a positive integer, (through 
the set of QP peaks defined by Eq. (|12|l ). to another set 
of periodic peaks <I>'j = m^[, when the "quasiperiodic- 
ity parameter" 7 = as/aL is gradually tuned between 
the values 7 = 1.0 and 7 = (see Fig. 2). These limits 
(7 = 1.0 and 7 = 0) correspond to periodic chains. 

Fig. 5 illustrates the corresponding evolution of the 
Fourier-transform of the distribution function of vortices 
pinned on a ID QP array of pinning sites. 

For 7 = 1 (Fig. 5a), there is a single sharp peak 
(accompanied by negligibly small satellites) correspond- 
ing to a periodic chain. For smaller 7's (e.g., 7 = 0.8 
(Fig. 5b), 7 = 0.8 (Fig. 5c), or 7 = 0.618 (Fig. 4a)), a 
set of satellite peaks appears around the main peak. Si- 
multaneously, the intensity of this peak decreases giving 
rise to another main peak for smaller value of q (Fig. 5d), 
which corresponds to a periodic chain with a larger pe- 
riod. 



V. PINNING OF VORTICES BY 2D 
QUASIPERIODIC PINNING ARRAYS 

In the previous section we studied the pinning of vor- 
tices by ID QP chains of pinning sites. In particular, we 
showed how the quasiperiodicity manifests itself in the 
critical depinning current, Jc Ny^"^ /Ny, when increas- 
ing the applied magnetic flux, ^ Ny. We found that 
the positions of the peaks of the function Jc{^) are gov- 
erned by "harmonics" of long and short periods of the 
QP chain of pinning sites. Independently of the length 
of the chain (for Np > 21), the peaks form a QP sequence 
including the Fibonacci sequence as a fundamental sub- 
set. This self-similarity effect is clearly displayed in the 
Fourier-transform of the distribution function of the vor- 
tices on a ID QP array of pinning centers. The evolution 
of QP peaks, when gradually changing the "quasiperiod- 
icity" parameter 7, has revealed a continuous transition 
from a QP chain — through the set of QP states (most 
pronounced for 7 = 1/r « 0.618) — to another periodic 
chain, 7 = 0, with a longer period. This phenomenon has 
been studied both in real space and in reciprocal fc-space. 

In the present section and in the next sections, we an- 
alyze vortex pinning by 2D QP arrays; in particular, by 
an array of pinning sites placed in the nodes of a five- 
fold Penrose lattice. Before tackling the Penrose-lattice 
pinning array, let us start with a simplified system which 
one can call "2D-quasiperiodic" (2DQP) since it is a 2D 
system periodic in one direction (x-direction) and QP 
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in the other direction (y-direction). 

A. 2D-quasiperiodic triangular lattice of pinning 
sites 

In an infinitely long one-dimensional homogeneous su- 
perconductor without any pinning centers, vortices ob- 
viously are equidistantly distributed, forming a peri- 
odic chain. Similarly, as it is well-known, in a three- 
dimensional superconductor (or in quasi-two-dimensional 
slabs or films), vortices organize themselves in a peri- 
odic triangular lattice shown schematically in Fig. 6a. 
If we keep the lattice undistorted along the a;-direction 
and introduce a quasiperiodic deformation along the y- 
direction, similarly to the case of a ID QP chain, we 
obtain a 2QP triangular lattice as shown in Figs. 6b, c. 
The "quasiperiodicity" parameter 'jy for this 2DQP lat- 
tice is defined as the ratio of the short to long peri- 
ods, OL to as (see Fig. 6c). The 2DQP triangular ar- 
rays of pinning sites arc shown for 7j, = 1.0 (Fig. 6a), 
jy = 1/r, where t = {1 + V5)/2 (Fig. 6b), jy = 0.5 
(Fig. 6c). The corresponding functions Jc{Ny ^ $) are 
presented in Figs. 6d,e,f for the following pinning param- 
eters: /p//o = 2.0, and Vp = 0.1 A. 

For the triangular array of pinning sites (Fig. 6a), 
we obtain a well-resolved main commensurability peak 
(Fig. 6d), corresponding to the first matching field, at 
$ = $1. Note that the vortex lattice is in general in- 
commensurate with a triangular lattice of pinning sites 
for the second matching field, i.e. when <& = 2<I>i (see, 
e.g., Ref. "2^). For instance, for the parameters used in 
our simulations, only each second row of pinning sites is 
occupied at $ = 2<I>i, resulting in a very weak maximum 
of the function Jc(iV„ ~ $) at that point; the param- 
eters used {fp/fo = 2.0, rp ~ O.IA) are nearly optimal 
for revealing features of the function Jc{^) related to 
quasiperiodicity. 

When tuning 7 out of the periodic value jy = 1.0, 
the main peak decreases in magnitude, and a maximum 
forms near it at a larger value of A^t, ~ $. These changes 
are demonstrated in Fig. 6e (the corresponding pinning 
array is shown in Fig. 6b) for 7^ = 1/r. 

It should be noted that here the parameter "fy has a 
different meaning, in the case of triangular 2DQP lat- 
tices, compared to the case of 7 for the ID QP chains 
considered in the previous section. In a ID QP chain, 
7 — as/oL is the ratio of distances between pinning sites, 
whereas in a triangular 2DQP lattice 7y defines the ratio 
of the distances between the rows of pinning sites (see 
Fig. 6c) . It is easy to show that the ratio of distances be- 
tween the neighboring pinning sites in a triangular 2DQP 
lattice is 

where a is the period of the 2DQP triangular lattice along 



the (periodic) a:-direction. Thus, for jy — l/r, the pa- 
rameter 7' defined by Eq. (|14|l becomes 7' ~ 0.7. 

For 7j, = 0.5 (7' ~ 0.6), the function Jc(iV„) is plotted 
in Fig. 6f. The main peak is further depressed, while 
the closest satellite peak becomes more pronounced. In 
addition, other satellite peaks appear, which are much 
less pronounced. 



VI. TRANSITION FROM A TRIANGULAR TO 
A QUASIPERIODIC PENROSE-LATTICE 
ARRAY OF PINNING SITES 

Above, we have revealed some features of the behavior 
of the critical depinning current Jc as a function of the 
applied magnetic flux $ (or as a function of the number 
of vortices in the system ~ $). They have been 
found under the transformation of a triangular lattice to 
a 2DQP triangular lattice (array) of pinning sites. 

Consider now a similar procedure but the final configu- 
ration of the transformation from a triangular lattice will 
be a 2D QP array of pinning sites, namely, an array of 
pinning sites located at the nodes of a five-fold Penrose 
lattice. This kind of lattice is representative of a class 
of 2D QP structures, or quasicrystals, which are referred 
to as Penrose tilings. These structures possess a local 
order and a rotational (five- or ten-fold) symmetry, but 
do not have translational long-range order. Being con- 
structed of a series of building blocks of certain simple 
shapes combined according to specific local rules, these 
structures can extend to infinity without any defects |llj | . 
Below, we will discuss in more detail the structure of the 
Penrose lattice. 

The transformation of a triangular lattice to the Pen- 
rose one is a rather non-trivial procedure, as distinct 
from the transformation to a 2DQP triangular lattice 
done above when we simply stretched some of the inter- 
row distances and squeezed other ones according to the 
Fibonacci rules (Eq. (jHJ) for a one-dimensional QP lat- 
tice. In order to find intermediate configurations between 
the triangular lattice and the Penrose lattice, we employ 
the following approach. First, we place non- interacting 
vortices at the positions coinciding with the nodes of 
the Penrose lattice (these can be considered as pinning 
sites for vortices, which right afterwards are "switched 
off"). Then we let the vortices freely relax undergoing 
the vortex-vortex interaction force at low temperatures 
(and no pinning force). The vortices relax to their ground 
state, which is a triangular lattice. During the relax- 
ation process, we do a series of "snapshots" , recording 
the coordinates of the vortex configurations at different 
times. The sets of coordinates obtained are then used 
as coordinates of pinning sites. We arrange these sets 
in antichronological order to model a continuous transi- 
tion of the pinning array from its initial configuration, a 
triangular lattice, to its final configuration, a Penrose 
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FIG. 6: The spatial distribution of pinning sites for a triangular lattice (i.e., periodic) (a), and 2D-quasiperiodic (2DQP) 
triangular, i.e. periodic along one direction (the a:-direction) and QP along the other one (the y-direction) , (b) and (c). The 
parameter 7^ is defined as a ratio of short (as) to long (az,) periods in the y-direction, as shown in (c). The values of the 
parameter 7^ are: 7^, = 1.0 (a), 7^ = 1/t, where t = {1 + y/E)/2 (golden mean) (b), 7^ = 0.5 (c). The critical depinning current 
Jc, as a function of the number of vortices, ~ $, for triangular (d) and the 2DQP triangular lattices [shown in (b) and (c)] 
(e) and (f), correspondingly, for fp/fo = 2.0, = O.IA. 
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FIG. 7: (This figure is available in "png" format: "Penrose 
Fig 7. png"; color) Left column: Transformation of a trian- 
gular lattice of pinning sites to a five-fold Penrose lattice. 
The distributions of the pinning sites in the triangular lattice 
(shown by blue solid circles in (a) and by blue open circles 
in (b), (c), and (d), for comparison). Intermediate configu- 
rations (shown by red solid circles in (b) and (c)) between 
the triangular (a) and the Penrose lattice (shown by red solid 
circles in (d)). Right column: The corresponding critical de- 
pinning currents, Jc, as a function of the applied magnetic 
flux, $, for the pinning arrays shown in (a) to (d), respec- 
tively: for the triangular lattice (shown by blue solid line in 
(e) and by blue dashed lines in (f), (g), and (h), for com- 
parison); for the intermediate conflgurations (shown by red 
solid lines in (f) and (g)); for the Penrose lattice (shown by 
the red solid line in (h)). For the Penrose lattice case in (b), 
the drop in Jc{^) is an artifact of the boundary conditions. 
Namely, the Penrose lattices of pinning sites did not flt the 
square cell used in the simulations. Thus, the freely-moving 
vortices near the edges significantly decreased the value of Jc, 
especially near <I>i. This problem will be dealt separately in 
Fig. 14 and in Eq. l(T(ijl . 



lattice. 

In Figs. 7a to 7d, four of these configurations of pin- 
ning sites are shown. The triangular lattice is presented 
in Fig. 7a. Two intermediate configurations are shown 
in Figs. 7b, c. The pinning sites plotted in Fig. 7d are 
located on the vertices of a Penrose lattice. For compar- 
ison, a pinning array in the form of a triangular lattice 
(Fig. 7a) is also presented in Figs. 7b,c,d as open blue 
circles. The functions Jc($) calculated for the pinning 
arrays shown in Figs. 7a to d, are plotted, respectively, 
in Figs. 7e to h. The function Jc(<&) for the triangu- 
lar lattice (Fig. 7e) is also plotted, for comparison, in 
Figs. 7g,f,h as a blue dashed curve. 

The main commensurability peak related to the first 
matching field in a triangular lattice of pinning sites, ob- 
served at $ = $1 (Fig. 7e), turns out to be rather sta- 
ble with respect to moderate deformations of the lattice 
(Fig. 7f, see also Fig. 7b). It still has a maximum height 
in Fig. 7f, although it broadens. However, the depths of 
the valleys near the peak decreases by about 20 to 30 per 
cent. Two sharp peaks near $ = ^i/3 and — <f>i/6 
(Fig. 7e), related to the commensurability of the long- 
range order in a triangular lattice, disappear. With fur- 
ther deformation, e.g., for the pinning array shown in 
Fig. 7c, the main peak still remains but only about 80 
per cent of the vortices are pinned in this case. The func- 
tion Jc(^) becomes somewhat smoother, and it does not 
display any pronounced features (for $ > 2$i/3) except 
the main maximum. 

The transition to a Penrose-lattice array of pinning 
sites (Fig. 7d) is accompanied by the appearance of a 
specific fine structure of the function Jc($). Namely, 
two well-resolved features on the broad main maximum 
(Fig. 7h) are the most pronounced ones. Other, less 



pronounced, features will be discussed below for larger 
Penrose-lattice arrays. For large arrays, the function 
Jc(<i>) is much less affected by fluctuations related to 
the entrance of each single vortex in the system, which 
are significant for the small-size array shown in Fig. 7d 
{Np — 56). This small-size array is used here just as an 
illustration, for studying the transition form a periodic 
(triangular) to a QP (Penrose lattice) pinning site array. 
However, studying even a relatively small piece of a QP 
structure provides some useful information about prop- 
erties of the whole system based on the self-similarity of 
the lattice, which was revealed for ID QP chains in the 
previous section, and which will be demonstrated below 
for 2D QP structures. 

The Penrose and the 2DQP triangular lattice (Fig. 6f) 
both have an important similar feature: their Jc{^) has 
two nearby maxima. Thus, our previous analysis based 
on several alternative ways to continuously deform a pe- 
riodic lattice to a QP one shows that the features shown 
are hallmarks of QP pinning arrays. 

In the next section, the origin of the features observed 
will be explained on the basis of a detailed analysis of the 
structure and of the building blocks forming a five-fold 
Penrose lattice. Other, less pronounced, features will also 
be discussed. Some of them will be found in larger arrays 
in our numerical simulations. 



VII. ANALYSIS OF THE FINE STRUCTURES 

OF THE FUNCTION Jc($) IN A 
QUASIPERIODIC PENROSE-LATTICE ARRAY 
OF PINNING SITES 

The structure of a five-fold Penrose lattice is shown in 
Fig. 8. As an illustration, a five-fold symmetric fragment 
which consists of 46 points (nodes of the lattice) is pre- 
sented (Fig. 8a). According to specific rules, the points 
are connected by lines in order to display the structure of 
the Penrose lattice (compare, e.g., to Fig. 7d). The ele- 
mental building blocks are rhombuses with equal sides 
a snd angles which are multiples of 9 = 36°. There 
are rhombuses of two kinds forming the Penrose lattice 
(Fig. 8b): (i) those having angles 26 and W (so called 
"thick"; they are empty in Fig. 8), and (ii) rhombuses 
with angles and 46 (so called "thin" ; they are colored 
in orange in Fig. 8). 

Let us analyze the structure of the Penrose lattice from 
the point of view ot its pinning properties, when pinning 
sites are placed in the vertices of the lattice. In par- 
ticular, we are interested whether any specific matching 
effects can exist in this system between the pinning lat- 
tice and the interacting vortices, which define the critical 
depinning current at different values of the applied mag- 
netic field (i.e., the function Jc (<&)). 

On the one hand, QP (quasicrystalline) patterns are 
intrinsically incommensurate with the fiux lattice for any 
value of the magnetic field therefore, in 

contrast to periodic (e.g., triangular or square) pinning 
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FIG. 8: (Color) The structure of a five- fold Penrose lattice 
(a). The elemental building blocks are rhombuses with equal 
sides a and angles which are multiples oi 6 = 36°. There are 
rhombuses of two kinds: those having angles 26 and W (so 
called "thick"), and rhombuses with angles 6 and AO (so called 
"thin") (b). The distances between the nodes, i.e. the lengths 
of the diagonals of the rhombuses are: 1.176a (the short diag- 
onal of a thick rhombus); (l-|-\/5)a/2 = ra ~ 1.618a (the long 
diagonal of a thick rhombus); [(1 -I- \/5)/2 — l]a = (r — l)a ~ 
0.618a, (the short diagonal of a thin rhombus); 1.902a (the 
short diagonal of a thin rhombus) (b). 



arrays, one might a priori assume a lack of sharp peaks 
in Jc{^) for QP arrays of pinning sites. 

On the other hand, the existence of many periods in 
the Penrose lattice can lead to a hierarchy of matching 
effects for certain values of the applied magnetic field, 
resulting in strikingly-broad shapes for Jc($). 

In order to match the vortex lattice on an entire QP 
pinning array, the specific geometry of the elements which 
form the QP lattice is important as well as their arrange- 
ment, as distinct from the flux quantization effects and 
superconductor-to-normal phase boundaries for which 
the areas of the elements only plays a role 0, 0, 0| . 
As mentioned above, a five-fold Penrose lattice is con- 
structed of building blocks, or rhombuses, of two kinds. 
While the sides of the rhombuses are equal (denoted by 
a), the distances between the nodes (where we place pin- 



ning sites) are not equal (which is problematic for vor- 
tices). The lengths of the diagonals of the rhombuses 
are as follows (Fig. 8b): 1.176a (the short diagonal of a 
thick rhombus); (1 -I- \/5)a/2 = ra « 1.618a, where r is 
the golden mean (the long diagonal of a thick rhombus); 
[(l + V5)/2-l]a (T-l)a « 0.618a, (the short diagonal 
of a thin rhombus); 1.902a (the short diagonal of a thin 
rhombus). 

Based on this hierarchy of distances, we can predict 
matching effects (and corresponding features of the func- 
tion Jc{'^)) for the Penrose-lattice pinning array. 

First, we can expect that there is a "first matching 
field" (let us denote the corresponding flux as $i) when 
each pinning site is occupied by a vortex. Although sides 
of all the rhombuses are equal to each other similarly to 
that in a periodic lattice, nevertheless this matching ef- 
fect is not expected to be accompanied by a sharp peak. 
Instead, it is a broad maximum since it involves three 
kinds of local "commensurability" effects of the flux lat- 
tice: with the rhombus side a; with the short diagonal 
of a thick rhombus, 1.176a, which is close to a; and with 
the short diagonal of a thin rhombus, which is the golden 
mean times a, 0.618a (see Fig. 8b). 

It should be noted that this kind of matching assumes 
that a vortex lattice is rather weak, i.e. the effect can be 
more or less pronounced depending on the specific rela- 
tions between the vortex- vortex interaction constant and 
the strength of the pinning sites as well as on the dis- 
tance between pinning sites and their radius. Assuming 
that the vortex-vortex interaction constant is a material 
parameter, all others can be adjustable parameters in 
experiments with artificially-created QP pinning arrays. 
For instance, the pinning parameters can be "adjusted" 
by using as pinning centers antidots, i.e. microholes of 
different radii "drilled" in a superconductor film 
or blind antidots of different depths and radii. 

Further, we can deduce that next to the above "main" 
matching flux there is another matching related with the 
filling of all the pinning sites in the vertices of thick rhom- 
buses and only three out of four of the pinning sites in the 
vertices of thin rhombuses, i.e., one of the pinning sites in 
the vertices of thin rhombuses is empty. For this value of 
the flux, which is lower than $1, matching conditions are 
fulfilled for two close distances, a (the side of a rhombus) 
and 1.176a (the short diagonal of a thick rhombus) but 
are not fulfilled for the short diagonal, a/r, of the thin 
rhombus. 

Therefore, this QP feature is related to the golden 
mean value, although not in such a direct way as in the 
case of a ID QP pinning array. This 2D QP matching re- 
sults in a very wide maximum of the function Jc{^)- The 
position of this broad maximum, i.e., the specific value 
of $ (denoted here by $vacancy/thin = *v/t = 0.757 $1) 
could be found as follows. The ratio of the numbers 
of thick and thin rhombuses is determined by the Fi- 
bonacci numbers and in the limit of large pinning arrays, 
Np — > 00 this ratio tends to the golden mean. The num- 
ber of unoccupied pinning sites is governed by the num- 
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FIG. 9: (This figure is available in "png" format: "Pen- 
rose Fig 9. png"; color) The critical depinning current Jc as a 
function of the applied magnetic flux, $ ~ A'„, for an array 
of pinning sites placed at the nodes of a five-fold Penrose lat- 
tice (for a part of the lattice which contains A'j, = 46 pinning 
sites) (a) . The distributions of vortices (shown by green dots) 
pinned on the Penrose- lattice pinning site array (pinning sites 
are shown by red open circles connected by orange solid lines 
used in order to show the Penrose lattice structure, i.e., thick 
and thin rhombuses), for specific values of the applied mag- 
netic flux: (b) $ = $vacancy/thin = $v/t = 0.757 $i, VOrticCS 

occupy all the pinning sites except those in one of the two ver- 
tices (connected by the short diagonal) of each thin rhombus, 
each single and each pair of double thin rhombuses contain 
one unoccupied pinning site at the matching field $v/t; (c) 
$ = $1, the number of vortices 7V„ coincides with the num- 
ber of pinning sites A'j, , and almost all the vortices are pinned 
(because of using a square simulation cell, some of the vortices 
are always "interstitial" but allow to keep the average vortex 
density in the entire simulation cell; due to these additional 
vortices, the value of the function Jc($) effectively reduces by 
a "filling" factor -q = Ap/A « 0.575, where Ap and A are the 
areas of the Penrose-lattice "sample" and of the simulation 
region; (d) $ = "I>interstitiai/thick = '^i/T = 1.482 $1, vortices 
occupy all the pinning sites and interstitial positions inside 
each thick rhombus, one vortex per each thick rhombus. The 
parameters are /p//o ~ 2.0, Vp — O.IA. 

ber of thin rhombuses. However, some of the thin rhom- 
buses are separated from other thin rhombuses by thick 
ones (call them single thin rhombuses), but some of them 
have common sides with each other (double thin rhom- 
buses). Therefore, the number of vacancies (i.e., unoccu- 
pied pins) is then the number of single thin rhombuses 
plus one half of the number of "double" thin rhombuses, 

iV7(<i>v/t) = n:,, + i<, (15) 

where iV™ is the number of unoccupied pinning sites at 

$ = $v/t, -^rh ^^-d N^Y^ arc the numbers of single and 
double thin rhombuses, correspondingly. 

For higher vortex densities (e.g., for <i> = 
^'intorstitiai/thick = *i/T = 1.482 we Can expect 
the appearance of a feature (maximum) of the func- 
tion Jc(^) related to the entry of a single interstitial 
vortex into each thick rhombus. The position of this 
maximum is determined by the number of vortices at 
$ = $1, which is Nv{^) = Np, plus the number of 
thick rhombuses, N^^^^^^. Since the ratio of the number 
of thick to that of thin rhombuses is the golden mean, 
T ~ (1 + \/5)/2 (in an infinite lattice; in a finite lattice, 
it is determined by a ratio of two successive Fibonacci 
numbers), then N^^^'^^ = N^-h/r, where A^rh is the total 
number of rhombuses. Here we used: 1/r = r — 1. 

In Fig. 9a, the function Jd^) is plotted for an array of 
pinning sites in the form of a part of the Penrose lattice, 
shown in Figs. 9c,d,e, which consists of 20 thick {N^^^^^ = 



20) and 15 thin (Af*h'" = 15) rhombuses and contains 46 
nodes (pinning sites). The nodes are connected by lines 
in order to show the rhombuses. 

The distribution of vortices for $ = $1 is shown in 
Fig. 9c. The number of vortices Ny coincides with the 
number of pinning sites Np, and almost all the vortices 
are pinned. Note that since we use a square simulation 
cell, some of the vortices are always outside the "Penrose 
sample" . These vortices mimic the externally applied 
magnetic field and determine the average vortex density 
in the entire simulation cell. Because of these additional 
vortices, the value of the function Jc($) effectively re- 
duces approximately by a "filling" factor 77 which is 

Ap 

= -i- w 0.575. (16) 

Here Ap and A are the areas of the Penrose lattice (i.e., 
the area of all the rhombuses) and of the simulation re- 
gion. 

The value of the function Jc($) in the maximum 
$ = $1 (Fig. 9a) is Jc ~ 0.55, i.e. corresponds to al- 
most perfect matching (two pinning sites occurred to be 
unoccupied in the distribution shown in Fig. 9c) taking 
into account Eq. (|16|l . 

Let us now more carefully analyze the calculations of 
Jc{^) for the Penrose lattice, presented in Fig. 9. In 
Fig. 9b, the distribution of vortices is shown for $ = 
$v/t- Vortices occupy all the pinning sites except those 
situated in one of the two vertices, connected by the short 
diagonal, of each thin rhombus. Thus, each single and 
each pair of double thin rhombuses contain one vacancy 
(unoccupied pinning site) at the matching field $v/t- The 
corresponding maximum is indicated by the arrow (b) in 
Fig. 9a. 

The location of vortices for $ = $1 (the maximum 
(c) in Fig. 9a) is shown in Fig. 9c; here the number of 
vortices Ny coincides with the number of pinning sites 
Np, and almost all the vortices are pinned. 

The distribution of vortices for $ = *i>i/T (Fig- 9d) is 
also in agreement with our expectation: vortices occupy 
all the pinning sites (there is only a single "defect" in 
the distribution shown in Fig. 9d: one vortex left the 
pinning site and became interstitial) plus interstitial po- 
sitions inside each thick rhombus, i.e., one vortex per 
each thick rhombus. However, the corresponding feature 
of the function Jc{^) (arrow (d) in Fig. 9a) is less pro- 
nounced than the two above maxima at $ = $1 and at 

In addition, there is a weak feature of the function 
Jc{^) at $ w <I>v/t/2, which more clearly manifests itself 
for larger Penrose-lattice pinning arrays (see Fig. 10a). 

Therefore, the calculated distributions of the vortices 
pinned on the Penrose-lattice pinning site array and the 
resulting function Jc(<&) have revealed the QP features 
which are in agreement with our expectations. The spe- 
cific structure of the function Jd^) is consistent with 
two previous derivations both based on continuously de- 
forming a QP lattice into a Penrose one (Sections V and 
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FIG. 10: (This figure is available in "png" format: "Penrose 
Fig 10. png"; color) The critical depinning current Jc as a 
function of the applied magnetic flux, $ ~ A'^ , for an array of 
pinning sites placed in the nodes of a flve-fold Penrose lattice 
(for TVp — 301 ) (a). The distributions of vortices (shown 
by green solid circles) pinned on the Penrose-lattice pinning 
site array (shown by red open circles), for speciflc values of 
the applied magnetic flux which correspond to two matching 
fields: (b) $ = $v/t, vortices occupy all the pinning sites 
except one in each thin rhombus; (c) $ = $1, all the vortices 
are pinned. The parameters are the same as in Fig. 9. 



VI). 

In Fig. 10a, the function Jc{Nv) is shown for a larger 
Penrose-lattice array of pinning sites, Np = 301. The 
above QP features in Jc{Ny) are much more pronounced 
in this case then for smaller arrays because of a consid- 
erable reduction of the "noise" related with an entry of 
each single vortex in the system. 

In particular, the main maximum of the function 
JciNy), which corresponds to the matching condition 
$ = $1 {Ny = 301), transforms into a rather sharp peak 
with the magnitude ry. Also, a local maximum of Jc{^) 
at $ R:! $v/t/2, is more pronounced for Np — 301, as 
mentioned above. 

Finally, Fig. 11 demonstrates the function Jc{^), cal- 
culated for different samples with Np = 46, 141, and 
301 (Fig. 11a), and also for different criteria of Jc- for 
the "static" and dynamical criteria (Fig. lib). In the dy- 
namical simulations of Jc using a threshold criterion, i.e., 
Jc is obtained as the minimum current J (x ff which 
depins the vortices. In the Appendix A, we show on- 
set of vortex motion when the applied current J exceeds 
the critical currnet: J > Jc- The results obtained using 
these two criteria are essentially equal, and throughout 
this work we use the "static" criterion defined above. 



VIII. ANALYTICAL APPROACH 

The following analysis could provide a better under- 
standing of the above structure of Jc{^) for the Penrose 
pinning lattice. Let us compare the elastic Ed and pin- 
ning i^pin energies of the vortex lattice at Hi and at (the 
lower field) H^/t, corresponding to the two maxima of Jc 
(e.g.. Figs. 9a, 10a, 11). Vortices can be pinned if the 
^pin P "pin of the pinnin g energy i s larger 
than the increase of the elastic energy related 



gain Epin t/pin/3"pin 

than the increase of tl 
to local compressions 



E, 



Ci 



(17) 



The shear elastic energy (oc Cge) provides the same qual- 



itative result. Here, 
pinning centers. 



pm 
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is the density of 
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FIG. 11: (Color) (a) The critical depinning current Jc as a 

function of the applied magnetic fiux, $ ~ Ni^-' /Nv, shown 
for three different Penrose-lattice pinning site arrays: A'p = 
46 (shown by blue solid line), A'p = 141 (shown by green 
solid line), A'p = 301 (shown by red solid line). The revealed 
features are hallmarks of a Penrose-lattice pinning site array, 
(b) Comparison of Jc as a function of the applied magnetic 
fiux, calculated using the "static" and dynamical criteria. The 
latter means calculating Jc using a threshold criterion, i.e., Jc 
is obtained as the minimum current J oc }f which depins 
the vortices. The results obtained using these two criteria are 
essentially equal. 



and P{H > Hi) — 1 \s the fraction of occupied pinning 
sites (/? = 1 for iJ = and /3 = 0.757 for H = iJ^/t), 

1 /2 

floq = (2/-\/3/3npin) is the equilibrium distance be- 
tween vortices in the triangular lattice, b is the minimum 
distance between vortices in the distorted pinned vortex 



lattice (6 
and 



a/r for H = Hi and 5 = a for iJ = H^u), 



Cii = 



47r(H-A2fc2) 



(18) 



PiH < iJi) = H/Hi = S/($onpi„), 



is the compressibility modulus for short-range deforma- 
tions |3l| with characteristic spatial scale k « (jT-pin)^^^- 
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The dimensionless difference of the pinning and elastic 
energies is 



47rA2 



where 



J diff — ^2 " 



$2 



1-6 



/3\/3npi, 



1/2- 



(19) 



(20) 



The function /diffO is shown schematicaUy in the inset 
to Fig. 12. Near matching fields, Jc has a peak when 
/diff > (and no peak when fdis < 0). Since only two 
matching fields provide /dig > 0, then our analysis ex- 
plains the two-peak structure observed in Jc shown in 
Figs. 9a, 10a, 11. For instance, for the main matching 
fields Eq. m gives: /diff($v/t) ~ 0.0056, /diff($i) « 
0.0058, and /diff($i/T) ~ -0.09. Note that for weaker 
pinning, the two-peak structure gradually turns into one 
very broad peak, and eventually zero peaks for weak 
enough pinning (see Fig. 12). The Jc peaks correspond- 
ing to higher matching fields are strongly suppressed be- 
cause of the fast increase (oc B^) of the compressibil- 
ity modulus Cii and, thus, the elastic energy with re- 
spect to the pinning energy; the latter cannot exceed the 



maximum value Uu 



The subharmonic peaks of Jc , 



which could occur for lower fields H < H, 



v/t! 



are also 



suppressed due to the increase of Cn associated with the 
growing spatial scales 1/fc of the deformations. 



IX. THE CRITICAL CURRENT Jc($) IN A 
RANDOMLY DISTORTED TRIANGULAR 
LATTICE 

Above we have studied the function Jc($) for periodic, 
QP 2D arrays of pinning sites and analyzed the transition 
from the periodic triangular lattice to the QP Penrose 
lattice (see Fig. 7). One of the issues, which is related 
to this analysis and can be useful for practical applica- 
tions, is the increase of the critical current (shown, e.g., 
in Fig. 7f) in the regions corresponding to minima of 
Jc($) for periodic (triangular) pinning arrays. The situ- 
ation shown in Fig. 7f seems to be the optimal from the 
point of view of a homogeneous increase of Jci'^) with- 
out degradation of the main peak at <i> = $i. Recall that 
it corresponds to a slightly "quasiperiodically distorted" 
triangular lattice (see Fig. 7b). The pinning sites of the 
triangular lattice are shifted from their "correct" posi- 
tions but not randomly: their positions are determined 
by vortex-vortex interaction, which tries to restore the 
triangular lattice, and by the memory about previous 
configurations including the initial one, i.e., the Penrose 
lattice. Analyzing the lattice presented in Fig. 7b we can 
deduce that it keeps some short-range order (i.e., dis- 
torted triangular cells similar in shape and size to those 
in the triangular lattice) but does not have long-range or- 
der of the triangular lattice. As a result, the main peak 
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FIG. 12: (Color) The critical current Jc("l>) for Penrose- 
lattice arrays for different pinning strength /p//o = 2.0 (red 
solid line), fp/fo — 1.4 (green solid line), fp/fo = 0.5 (blue 
solid line). The peak at "I>i is suppressed for weaker pinning 
(fp/fo = 1.4). Eventually, all the main peaks disappear for 
sufRciently weak pins {fp/fo ~ 0.5). The inset shows the 
dimensionless difference, /diff, of the pinning and the elastic 
energies versus the pinning-to-interaction energy ratio, for the 
broad Jc peak at "l>v/t (red dashed line) and for "l>i (red solid 
line). Only /diff > gives stable peaks in Jc- 



remains, since it is related to the short-order matching 
effects (i.e., over the distances of the order inter-site spac- 
ings a). The sharp decrease around the maximum and 
appearance of the deep valleys is explained by the ab- 
sence (due to long-range order, i.e., over distances longer 
than a) of any matching effects for the flux densities 
close to $ = $1 (there is no matching for $ — <i>i/2 or 
$ = 1.5$i and $ = 2$i for the triangular lattice). When 
the long-range order is destroyed, as shown in Fig. 7b, 
matching effects other than <I> = $i become allowed. 

It is appropriate to mention here that the QP Pen- 
rose lattice possesses a short-range order but does not 
have a long-range translational order. In such a way, the 
"quasiperiodically distorted" triangular lattice (Fig. 7b) 
or the QP Penrose lattice itself are good candidates for 
the optimal enhancement of the critical current in the 
regions where the function Jc{^) have minima for a pe- 
riodic lattice. (It should be noted here that the curves 
shown, e.g., in Fig. 7h, for the triangular and the Pen- 
rose lattices are calculated for the same cell, although 
the effective area of the Penrose lattice is smaller than 
that of the triangular lattice with the same number of 
pinning sites. This discrepancy is taken into account by 
the "filling factor" introduced by Eq. llt)|) . It should be 
also recalled, when comparing the function Jd^) for the 
case of the triangular and the Penrose lattices, that the 
main maximum of the curve for the Penrose lattice (see 
Fig. 11) is the second sharp peak at $ = $i.) 

In this respect, it is interesting to compare the above 
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FIG. 13: (This figure is available in "png" format: "Penrose 
Fig 13. png"; color) Left column: Randomly distorted triangu- 
lar lattices for r™^ = 0.2(a/2) (a), 0.3(a/2) (b), 0.4(a/2) (c), 
and for a/2 (d). For comparison, the triangular lattice is also 
shown. Right column: The functions Jc(*&) corresponding to 
the distributions of pinning sites shown in (a) to (d), are pre- 
sented, respectively, in (e) to (h). At low levels of noise (e), 
(f), the valleys (minima) start to fill due to disappearance of 
the long-range order, similarly to the case of the Penrose lat- 
tice, although accompanied with a weaker enhancement of Jc 
than for the case of the "quasiperiodic distortion" (Fig. 7f). 
For higher levels of noise, the main peak degrades without 
any essential enhancement of Jc in the neighborhood (g), (h). 




FIG. 14: (Color) The critical current Jc($) for a 301-sites 
Penrose-lattice (dark green solid line), (recalculated for flux 
only on the Penrose area, Ap), triangular (black solid line) 
and random (red solid line) pinning arrays. The Penrose lat- 
tice provides a remarkable enhancement of Jc($) over a very 
wide range of values of ^ because it contains many periods in 
it. 



results for the "quasiperiodic distortion" of the triangular 
lattice with its random distortion. For this purpose, we 
introduce a random angle aran: < aran < Stt, and a 
random radius of the displacement dran: < rj-an < t"'^, 
where r™^'^ is the maximal displacement radius, which is 
a measure of noise measured in units of a/2, where a is 
the (triangular) lattice constant. 

In Fig. 13, randomly distorted triangular lattices are 
shown for r™^'^ = 0.2(a/2) (Fig. 13a), 0.3(a/2) (Fig. 13b), 
0.4(a/2) (Fig. 13c), and for a/2 (Fig. 13d). For compar- 
ison, the triangular lattice is also shown in Figs. 14a to 
d. The corresponding functions Jc{^) are presented, re- 
spectively, in Figs. 13e to h. At low levels of noise (e.g., 
Fig. 13e, f) the valleys (minima) start to fill due to the 
disappearance of the long-range order, similarly to the 
case of the Penrose lattice, although accompanied with 



a weaker enhancement of Jc than for the case of the 
"quasiperiodic distortion" (Fig. 7f). For higher levels 
of noise, the main peak degrades without any essential 
enhancement of Jc in the neighborhood (Figs. 13g, h). 

For comparison, we also show in Fig. 14 the Jc{^) for a 
Penrose-lattice (calculated for the sample with Np = 301, 
only for the area of the Penrose lattice Ap, see Eq. H16|l ). 
Notice that the QP lattice leads to a very broad and po- 
tentially useful enhancement oi the critical current Jc{^), 
even compared to the triangular or random APS. The re- 
markably broad maximum in Jd^) is due to the fact that 
the Penrose lattice has many (infinite , in the thermody- 
namic limit) periodicities built in it |ll| . In principle, 
each one of these periods provides a peak in Jc(<i>). In 
practice, like in quasicrystalline difraction patterns, only 
few peaks sre strong. This is also consistent with our 
study. Furthermore, the pinning parameters can be ad- 
justed by using as pinning centers either antidots (mi- 
croholes of different radii "drilled" in the film P, 0] , or 
blind antidots of different depths and radii. Thus, our 
results could be observed experimentally. 



X. CONCLUSIONS 

The critical depinning current Jc, as a function of the 
applied magnetic flux $, has been studied in QP pinning 
arrays, from one-dimensional chains to two-dimensional 
arrays of pinning sites set in the nodes of quasiperiodic 
lattices including a 2D-quasiperiodic triangular lattice 
and a five-fold Penrose lattice. 

In a ID quasiperiodic chain of pinning sites, posi- 
tions of the peaks of the function Jc($) are governed by 
"harmonics" of long and short periods of the quasiperi- 
odic chain. Independently of the length of the chain, 
the peaks form a set of quasiperiodic sequence includ- 
ing a Fibonacci sequence as a basic subset. Analyzing 
the evolution of the peaks, when a continuous transi- 
tion is performed from a periodic to a quasiperiodic lat- 
tice of the pinning sites, we found that the peaks re- 
lated to the Fibonacci sequence are most pronounced 
when the ratio of lengths of the long and the short pe- 
riods is the golden mean. A comparison of the sets 
of peaks for different chains shows that the functions 
Jc{^) for the ID quasiperiodic chain is self-similar. In 
the /c-space, the self-similarity effect is displayed in the 
Fourier-transform of the distribution function of the sys- 
tem of vortices pinned on a ID quasiperiodic array of pin- 
ning centers. The evolution of quasiperiodic peaks when 
gradually changing the "quasiperiodicity" parameter 7 
(i.e., ratio of the lengths of short to long elements of a 
quasiperiodic chain) has revealed a continuous transition 
from a periodic chain - through the set of quasiperiodic 
states - to another periodic chain with a longer period. 
This phenomena has been studied both in real space and 
in reciprocal fc-space. 

In 2D quasiperiodic pinning arrays (e.g., Penrose lat- 
tice), the pinning of vortices is related to matching 
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conditions between a triangular vortex lattice and the 
quasiperiodic lattice of the pinning centers. Although 
more complicated than in ID pinning chains, the specific 
behavior of Jc($) is determined by the presence of two 
different kinds of elements - thick and thin rhombuses - 
forming the quasiperiodic lattice. Based on these consid- 
erations, the positions of the main maxima of Jc(^) for 
Penrose lattice are predicted. 

In particular, for the first matching field each pinning 
site is occupied by a vortex. The corresponding maxi- 
mum of the function Jc($) is broad since it involves at 
least three kinds of local matching effects of the flux lat- 
tice, with the rhombus side and with short diagonals of 
thick and thin rhombuses. 

Another Penrose-lattice matching field is related with 
local matching effects which involve the intervortex dis- 
tance of the vortex flux lattice, the rhombus side and the 
short diagonal of thick rhombus. For this field, all the 
pinning sites are occupied, which are situated in the ver- 
tices of thick rhombuses and only three out of four in the 
vertices of thin rhombuses. The number of unoccupied 
pinning sites is governed by the number of thin rhom- 
buses. Some of the thin rhombuses are single (i.e., sep- 
arated from other thin rhombuses by thick ones), while 
some of them arc doiiblc (i.e., have common sides with 
each other). Therefore, the number of vacancies is the 
number of single thin rhombuses plus one half of the 
number of "double" thin rhombuses. One more impor- 
tant feature of the function Jc(^) occurs for higher vor- 
tex densities, when a single interstitial vortex enters each 
thick rhombus. 

Numerical simulations performed for various sample 
sizes have revealed a good agreement with our predic- 
tions. 

The revealed features can be more or less pronounced 
depending on specific relations between the vortex- vortex 
interaction constant and the strength of the pinning sites, 
as well as on the distance between pinning sites and their 
radius. While the vortex- vortex interaction constant is a 
material parameter, all others can be adjustable param- 
eters in experiments with artificially created quasiperi- 
odic pinning arrays. This can be reached by using, for 
instance, antidots (i.e., microholes "drilled" in a super- 
conductor film) or blind antidots of different depths and 
radii as pinning centers. Our calculations provide the 
necessary relations between these parameters for possi- 
ble experimental realizations. 

A continuous deformation of the Penrose lattice to a 
periodic triangular lattice (i) shows that the above re- 
vealed features are hallmarks of quasiperiodic pinning ar- 
rays; (ii) provides us with a tool for the controlled change 
of the magnitude, sharpness and the position of the peaks 
of Jc(^) that is important for possible applications. In 
particular, our analysis shows that the quasiperiodic lat- 
tice provides an unusually broad critical current Jc($), 



that could be useful for practical applications demanding 
high Jc's over a wide range of fields. 



ACKNOWLEDGMENT 

This work was supported in part by the National Se- 
curity Agency (NSA) and Advanced Research and De- 
velopment Activity (ARDA) under Air Force Office of 
Scientific Research (AFOSR) contract number F49620- 
02-1-0334; and also supported by the US National Sci- 
ence Foundation grant No. EIA-0130383, and RIKEN's 
President's funds. 



APPENDIX A: Onset of vortex motion for currents 
higher then the critical current: J > Jc 

Here we present vortex flow patterns for currents ex- 
ceeding the critical value, J > Jc- In Fig. 15, the vortex 
flow patterns are shown calculated for a Penrose sample 
with Np = 141 (pinning sites shown by red circles). The 
ground-state vortex configuration is shown in Fig. 15a 
(vortices shown by green dots) for $ ft: $i and when 
no driving force is applied. This vortex configuration is 
similar to those shown in Figs. 9c and 10c, when the num- 
ber of vortices (within the sample area) is equal to the 
number of pinning sites, and all the vortices arc pinned. 
When an increasing driving force /d ~ J is applied, the 
vortices do not move until fd reaches some threshold 
value, when vortices depin. The current which corre- 
sponds to the driving force depinning the vortices, is then 
defined as the critical current Jc (dynamical criterion). A 
comparison of Jc's calculated using this criterion and the 
static criterion, is shown above in Fig. lib. When un- 
pinned, the vortices move along complicated trajectories, 
or "channels" created by pinning arrays and interacting 
with other vortices. Figs. 15b to 15d show the onset of 
the flux motion for J > Jc, following the traces of mov- 
ing vortices over distances about 0.5A (b), lA (c), and 
2A (d). On the subsequent consecutive snapshots, vor- 
tices are shown by a sequence of conseciitivc open blue 
circles (and by blue solid circles for the last snapshot on 
each panel). These show dynamical configurations of the 
vortex lattice in motion. Note the appearance of local 
"rivers" of vortices moving along the channels between 
neighboring pinning sites. 
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FIG. 15: (This figure is available in "png" format: "Penrose 
Fig 15. png" ; color) Vortex flow patterns for J > Jc, calculated 
for a Penrose sample with Np = 141 (pinning sites shown by 
red circles), (a) Ground-state vortex configuration (vortices 
shown by green dots) when $ « $i, and no driving force 
is applied, (b-d) The onset of the flux motion for J > Jc, 
following the traces of moving vortices over distances about 
0.5A (b), lA (c), and 2A (d). On the subsequent consecutive 
snapshots, vortex trajectories are shown by black dotted lines 
(the blue solid circles show the last snapshot on each panel). 
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